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■ Abstract. Because we know little about the Galactic force-field away 

from the plane, the Galactic mass distribution is very ill-determined. 
I show that a microlensing survey of galaxies closer than 50Mpc would 
enable us to map in three dimensions the Galactic density of stellar mass, 
which should be strictly less than the total mass density. A lower limit can 
^ ! be placed on the stellar mass needed aX. R < Rq to generate the measured 

optical depth towards sources in the bulge. If the Galaxy is barred, this 
limit is lower by a factor of up to two than in the axisymmetric case. 
^> ' Even our limited knowledge of the Galactic force field suffices to rule 

CN ; out the presence of the amount of mass an axisymmetric Galaxy needs to 

' generate the measured optical depth. Several lines of argument imply that 

. the Galaxy is strongly barred only at R < A kpc, and if this is the case, 

' even barred Galaxy models cannot generate the measured optical depth 

. without violating some constraint on the Galactic force-field. Galactic 

I mass models that are based on the assumption that light traces mass, for 

. which there is significant support in the inner Galaxy, yield microlensing 

Qh| optical depths that are smaller than the measured value by a factor of 

Q . more than 2.5. 

> 

^ ' 1. Introduction 



Notwithstanding the difficulties to which heavy obscuration by dust gives rise, 
we now have a reasonable idea of how luminosity is distributed in the Galaxy. 
We know we live in a galaxy that has a bar, 3 — 4 kpc long, whose nearer 
end lies at positive longitudes. Within the bar there is a moderately flattened 
component in which the luminosity density rises roughly as r~^'^ with decreasing 
galactocentric radius r. Around the bar there is a roughly exponential disk with 
scale length ~ 3 kpc, and a moderately flattened metal-poor halo in which the 
luminosity density varies as ~ r~^-^. 

The situation regarding the Galactic distribution of mass is very much less 
satisfactory. Observations of external galaxies and cosmological theory have 
convinced us that light is not a good tracer of mass. Hence we cannot straight- 
forwardly translate our models of the luminosity distribution into models of the 
mass distribution; in principle, we should start afresh and derive the mass dis- 
tribution by tracing out the Galaxy's gravitational field F(r) and then applying 
the divergence operator: V • F = —A'KGp. 
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From observations of gas that flows in the Galaetie plane we have a fair 
idea of what F(r) looks like at points r within the plane, but we have very little 
secure knowledge of F out of the plane. This deficiency is serious, because until 
we know F away from the plane, we cannot apply V anywhere, and can say 
nothing with security about p anywhere] 

What we currently know about F out of the plane relates to points near the 
Sun. Creze et al. (1998) and Holmberg &; Flynn (2000) have used proper motions 
of Hipparcos stars to estimate the gradient of Fz near the plane and hence infer 
that the local mass density is 0.076 ± 0.015 pc"-'^ and 0.10 ± 0.01 M© pc'^, 
respectively, while Kuijken & Gilmore (1991) have used the radial velocities of 
stars at the SGP to estimate Fz at z = — l.lkpc and thus infer that the surface 
density of the Galaxy within 1.1 kpc of the plane is 71 zt 6 Mq pc~^. To constrain 
F at points that lie far from the Sun, we have to study objects that move to 
such points. High-velocity stars are the obvious tracers to use, because they are 
so numerous. They strongly constrain models of F because stars on essentially 
the same orbit can be studied both locally and in situ (Binney, 1994; Dehnen 
& Binney, 1996). Unfortunately, the potential of this approach has yet to be 
systematically exploited. Recently, there has been considerable interest in using 
tidal streamers associated with disrupted satellites as tracers of F (Johnston et 
al. 1999, Helmi et al. 1999). My own view is that high- velocity stars have greater 
potential because (a) they are vastly more numerous and (b) they do not require 
approximations of the level of the (manifestly false) assumption that all elements 
of a streamer are on the same orbit. Moreover, extracting useful information 
from streamers requires space-based astrometry, the requisite observations high- 
velocity stars are available now, and all that's lacking is machinery for modelling 
them. 

Gravitational microlensing directly probes the Galaxy's mass distribution, 
but in a very different way from classical studies of gas and stars. In fact, 
microlensing does not measure the smooth Galactic force-field F hut grainincss 
in the mass distribution. Hence, it is insensitive to the contribution to the 
latter from elementary particles and gas, and is therefore complementary to the 
traditional approach to the determination of p from F, rather than competitive 
with it. 

2. From optical depth to stellar density 

The optical depth to microlensing of a stellar object at distance sq is 



is the harmonic mean of the source-lens and lens-observer distances. If the source 
is extragalactic, s ~ s, and r becomes proportional to / ds p*s in the direction 
of the source. Suppose we measure r along many lines of sight all over the sky. 
Can we reconstruct p* from the data? 




(1) 



where 




(2) 
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Figure 1. Contours of equal microlensing optical depth for a distant 
source for two Galactic models. There are two contours per decade and 
the heavy contour is for r = 10'^. The upper panel is for the model 
one obtains from the luminosity model of Kent, Dame &: Fazio (1991) 
with an assumed mass-to-light ratio Tk = 1- The lower panel is for 
Model 1 of Dehnen &; Binney (1998), and all components, including the 
dark halo, have been assumed to contribute fully to r. Consequently, 
the optical depth at large |6| is larger in the lower than in the upper 
panel. 
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With some simplifying assumptions, we can. Binney &: Gerhard (1996) show 
that if the Galactic luminosity density j(r) is symmetric about the Galactic plane 
and two other, orthogonal planes (as it would be if it were triaxially ellipsoidal), 
then a Richardson-Lucy algorithm can be used to recover j(r) from its line- 
of-sight projection, = Jdsj, at each point 0, on the celestial sphere. It 

is straightforward to modify the derivation of Binney &: Gerhard to show that, 
with the same assumptions regarding symmetry, a Richardson-Lucy algorithm 
for the recovery of from T{i}) is 



P* 



(fc+i) 



«=^i"«i:^77^/i:^- (3) 



W(17,)s(r.)/ s(r,)' 



Here the sum over i is over the eight points rj that are connected to r by the 
assumed symmetry of the lines of sight to these points are in the directions 
fli, and their distances from the Sun are s(rj). 

While it is in principle possible to recover from r, it is not clear that 
this will ever be done. The problem is the small numerical factor in front of the 
integral in (1). Numerically, 



r = 6 X 10" 



IQio Mo/kpc 



(4) 



where we have again assumed that the source is extragalactic. Figure [H plots 
r for two typical Galactic models. The model underlying the upper panel is 
the iT-band luminosity model of Kent, Dame & Fazio (1991), which has been 
converted into a model of the stellar mass distribution by assuming a mass-to- 
light ratio Tk = 1, which appears to be correct for the solar neighbourhood 
(§10.4.4 of Binney & Merrifield 1998). One sees that r exceeds 10~^ only for 
lines of sight at fairly low latitudes. 

Existing data indicate that the duration of a microlensing event is typically 
tens of days, so in an observing season one obtains at most a handful of statis- 
tically independent observations per line of sight. Hence, in a given area of the 
sky the number of lines of sight that must be monitored for of order years to 
distinguish r from zero is ~ 1/r ^ 10^. Finding this number of extragalactic 
sources in each of a large number of patches of the sky is hard. Probably our 
best chance is offered by 'pixel lensing' towards nearby galaxies. Gould (1996) 
gives the necessary theory and I adopt his notation. A large role is played in this 
by the effective stellar flux , which is defined in terms of the stellar luminosity 
function (j) by 

F, = jdF<t){F)F^I jdF<t){F)F. (5) 

Per resolution element on a galaxy image, the rate at which microlensing events 
can be detected above a signal-to-noise threshold Qmin is 

^ - ^raF,, (6) 



Q 



mm 



where 1 > k,^ are dimensionless functions, and aF* is the rate at which the 
telescope detects photons from an object of flux F^. In terms of the surface 
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brightnesses of galaxy and sky S and 5sky, we have k = (l + 5sky/5')~"^ ~ S/Sgky 
over most of a galactic image. The value of depends on the degree to which 
the image is resolved into stars: if a significant part of the integral on the top 
of equation (P) comes from stars bright enough that lensing of them can be 
detected even at large impact parameter {u ^ 0.25), then ^ is small, with ^ ~ 1 
otherwise. 

Gould shows that if one is interested only in measuring r regardless of the 
masses of the lenses that generate it, the optimal observational strategy is to 
work in the regime ^ ~ 1 in which one detects only high-magnification events. 
We shall be in this regime provided 

where lv is the rate constant of a typical microlensing event and -Fpsf is the 
flux in a resolution element of the galactic image. How big will these quanti- 
ties be in a typical case? Suppose we are using a diffraction-limited telescope of 
diameter D to study a galaxy of distance s at radius R25 , where the l^-band sur- 
face brightness will be ~ 24magarcsec~'^. Gould gives the absolute magnitude 
corresponding to F^, as 

M^i = -4.84 + 3{V - I) (8) 

and estimates that from a star with / = 20 our telescope will collect photons at 
a rate 10(-D/1 m)^ s~^. Hence, 

aF, = 10(JLYl0-0Ai30-20-4.84+3iV-I))(^^Y\ (g^ 
Vim/ VlOMpc/ ^ ^ 

The I-band flux in the telescope's resolution element is 

Fpsf = FolO-°-<^^-(^-^)) [0.206(^/^^)1 ^ (10) 



F>/lm 

where Fq is some universal constant. We can express F* in terms of this same 
constant and M*/ thus 

VlOMpc^ ■ ^ ' 

When we substitute equations ( [lO[ ) and (|ll]) into equation and assume 
(Table 4.4 of Binney & Merrifield, 1998) 

K ^ 10-o.4(4-(V-7)) (12) 

and y — / ~ 1.25 (de Jong, 1995) the condition to be in the high-magnification 
regime becomes 

10.4f^)V^— )-V^^)-'(u;lweek)-i < (13) 
Vim/ VlOMpc/ VlOOOnm/ ^ ' An ^ ^ 
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Gould finds that Qmin ~ 7 is required for detection, so this condition is satisfied 
for ^ 4 m and s ~ 50Mpc. By equations (P), (^) and ( [l2|) the event rate in 
this regime is 



2^ . ^Q-0.4(6.66+2(y-/)) (J^\ V ^ ^ -2 g- 



Q^i, vim/ VlOMpc 

262r f D \2y s \ -2 



. f— ) f^— ) week-i (14) 
g^i^Vlm/ VlOMpc/ ^ ^ 

Hence, by directing a 4m telescope along a line of sight with r = 10^^, we should 
be able to detect three to four events per 10^ resolution elements per week. In 
fact a somewhat higher event rate could be achieved if the telescope were in 
space, because the event rate is inversely proportional to the sky brightness, 
which will at least 2 mag fainter in space than the ground-based value I have 
assumed. 

If the PSF of the telescope has FWHM of x arcsec, a galaxy of diameter 
y arcmin offers 3600(y/x)^ pixels to monitor. An L* galaxy at a distance of 
50Mpc has y ~ 1, and will provide 10^ pixels if x ~ 0.06. Thus, a diffraction- 
limited telescope of modest aperture monitoring galaxies within 50 Mpc for of 
order a year could map r sufficiently extensively for it to be possible to recover 
from a scheme such as (^) . 

Of course, the values of r recovered from such a survey would include contri- 
butions from self-lensing within the target object in addition to the optical depth 
through the Milky Way. As Crotts (1992) pointed out in relation to pixel lensing 
of M31, variations in r across the image of an highly inclined galaxy would help 
one to determine how much self lensing occurs in galactic halos rather than in 
disks or bulges. Consequently, the analysis of the data from the survey would 
proceed by modelling in some detail the distribution of r in each target galaxy, 
with the Milky Way's contribution in that direction as a single number to be 
fitted to a considerable body of data. 



3. Real data 

Before we lobby NASA and ESA for a dedicated microlensing space telescope, 
we should ask what can be learned from the existing microlensing data. Three 
areas of the sky have been monitored: (i) towards the Galactic bulge; (ii) towards 
some spiral arms; and (iii) towards the Magellanic Clouds. 

If one believes that the dark halo is comprised of elementary particles, only 
a very small optical depth is predicted towards the Clouds (~ 4 x 10~*). At 
1.2j^Q3 X 10"'' the measured optical depth is about three times larger, but still 
deriving from only 17 events (Alcock et al. 2000a). The nature of these events is 
controversial. There is a powerful case that many of the lenses lie in the Clouds 
themselves (Kerins &: Evans, 1999) because the only two lenses with reasonably 
securely determined distances (from finite-source effects) do lie in the Clouds. On 
the other hand, Ibata et al. (1999) may have detected a substantial population 
of old white dwarfs from their proper motions in the HDF. It is just possible 
that these objects provide the lenses for a significant number of the observed 
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events. The arguments against a large column density in white dwarfs remain 
powerful, however (Gibson & Mould 1997). 

The EROS collaboration (2000) has observed seven probable microlensing 
events in the directions of spiral arms, and from them derived an optical depth 
4.5lf^ X 10"'^ that is compatible with conventional Galaxy models. Unfortu- 
nately, the error bar on this measurement is rather large. 

3.1. The Galactic centre 

Several hundred events have been detected by various groups along lines of sight 
towards the Galactic centre, and these data pose a fascinating puzzle. They are 
harder to interpret than the data for lines of sight to the Clouds because we 
do not know a priori where the sources are. However, if we confine ourselves 
to the data for red-clump stars, we can have quite a precise estimate of their 
distribution down each line of sight. 

The red-clump stars must follow the general distribution of near infra-red 
light quite closely, because they are part of the population of evolved stars that 
are responsible for most of the Galaxy's near-IR luminosity. The DIRBE ex- 
periment aboard the COBE satellite mapped the Galaxy's IR surface brightness 
in several wavebands, and in the far-IR, where emission by dust is dominant. 
Spergel, Malhotra & Blitz (1996) used these data to estimate the effects of ex- 
tinction on the near-IR data, and produce maps of what the Galaxy would look 
like in the near-IR bands in the absence of extinction. 

Binney, Gerhard &: Spergel (1997) used their Richardson-Lucy algorithm 
to deproject these corrected near-IR data under the assumption that the Galaxy 
has three specified perpendicular planes of mirror symmetry. The upper panel 
of Figure |2| shows the luminosity density that they recovered when projected 
perpendicular to the plane from z = 225 pc upwards. Excluding the Galactic 
plane from the projection suppresses local maxima at ~ 3kpc along the y that 
are a notable feature of the density distribution in the plane (see Fig. 5 of Binney 
et al.). Binney et al argued that these maxima were artifacts resulting from the 
presence in the Galactic plane of spiral arms, which violate the assumed eight- 
fold symmetry. 

Bissantz & Gerhard (2000) have recently deprojected the same data with an 
entirely different technique. Rather than using the Richardson-Lucy algorithm, 
they formulate the deprojection problem as a regularized likelihood maximiza- 
tion. They do not directly impose any symmetry on the model but have a term 
in the penalty function that discourages deviations from eight-fold symmetry. 
Another term in the penalty function encourages luminosity to lie along the 
spiral arms delineated by Ortiz & Lepine (1993). 

The lower panel in Figure ^ shows that explicitly modelling the Galaxy's 
spiral structure in this way has the effect of making the bar longer and thinner 
than that recovered by Binney et al. (1997) - the axis ratio in the plane increases 
from 2:1 to 3:1. This change to the model bar enables the latter to reproduce an 
important datum that the Binney et al bar did not reproduce: the histograms 
from Stanek et al. (1994) that give for lines of sight at / ~ ±5° the number 
of clump stars at each apparent magnitude - see Figure ^. The ability of the 
model to reproduce these histograms to good accuracy strongly suggests that 
the model gives a faithful account of the distribution of red-clump stars. Hence, 
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Figure 2. Top panel: the Galaxy in the L band projected perpen- 
dicular to the plane from z = 225 pc upwards according to the model 
of Binney et al. (1997). Lower panel: the same view according to the 
model of Bissantz & Gerhard (2000). In both plots contours are log- 
arithmically spaced. In the upper plot there are three contours per 
decade, while below contours are explicitly labelled. 
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Figure 3. Predicted (curves) and measured (points) apparent- 
magnitude distribution of clump stars along three lines of sight through 
the bulge. The predictions are based on the model of Bissantz &; Ger- 
hard (2000). 



when this model is used to predict the microlensing optical depth to red-clump 
stars, we may be confident that any discrepancy does not arise from an incorrect 
distribution of source objects. 

Bissantz &; Gerhard convert their luminosity model of the inner Galaxy 
into a mass model by adopting the constant near-IR mass-to-light ratio T of 
Englmaier & Gerhard (1999). This value of T was obtained by comparing the 
pattern of gas flow predicted by the Binney et al. model for some assumed T 
to the observed {l,v) diagrams for HI and CO; if Bissantz Sz Gerhard were 
to repeat this exercise with their new photometric model, they would surely 
obtain a very similar value of T. Once T has been chosen, one can calculate 
the microlensing optical depth for red-clump stars along any line of sight. In 
Baade's window they obtain r = 1.24 x 10~^, which is essentially identical to the 
value obtained by Bissantz et al. (1997) from the shorter, fatter bar of Binney et 
al. (1997), and significantly short of the values from the MACHO collaboration: 
r = 3.91^ 2 X 10~^ (Id) directly measured for bulge clump giants by Alcock et al. 
(1997) and (3.9 ±0.6) x 10^^ estimated by Alcock et al, (2000b) for bulge stars 
from a difference-imaging analysis of an inhomogeneous collection of sources. 
Evidently one cannot obtain agreement with the MACHO optical depth under 
the assumption that mass follows light. 

Since the difference between the optical depth implied by constant T and the 
measured value is so large, it is natural to investigate an extreme model in which 
we ask, what is the minimum mass that is compatible with the red-clump optical 
depth attaining the MACHO value? Most of the red-clump sources lensed in 
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Baade's window {I, b) = (1°, —4°) lie close to the Galaxy's z axis, so we simplify 
the calculation by considering a source that lies distance h from plane on the 
axis. Consider the contribution to r from a band of mass M and radius r around 
the Galactic centre. If we assume that the band's surface density never increases 
with distance from the plane, then its mass will be minimized for a given optical 
depth when its surface density is constant and the line-of-sight to the source 
just cuts its edge. So we take the band's half-width to be h{RQ — r)/i?o, which 
makes the band's surface density 

Substituting this into equation (1) we find the band's optical depth to be 

_ GM 



(16) 



independent of radius (Kuijken, 1997). This minimum mass estimate holds if the 
mass is widely distributed in radius rather than concentrated in a single band, 
because we can imagine a radially continuous mass distribution to be made up 
of a large number of bands, and we have shown that when the band is optimally 
configured, its optical depth depends only on its mass. 

From equation (|l6|) (3.8 it 0.6) x lO^'' Mq is needed to produce the optical 
depth, to bulge sources that is implied by the latest MACHO results. A more 
realistic mass estimate is in excess of 8 x 10^*^ Mq because realistically we must 
assume that the surface density of the band falls off smoothly with distance 
from the plane, and if this decline is exponential with the optimal scale-height 
{h[l — r/i?o]), the band's mass must be e times that given by ( p!^ ) for a given 
optical depth, while if the vertical density profile is Gaussian with optimal scale- 



height {h[l — r/Ro]), equation (16) underestimates the band's mass by a factor 
\Afe72 ~ 2.07. 

For comparison, the mass of the Galaxy interior to the Sun is of order 
M = (220kms-^)2 x 8kpc/G ~ 8.9 x 10^° Mq. Thus this naivest estimate of 
the mass interior to the Sun is just barely equal to the minimum mass required 
in a circular configuration to produce the reported optical depth, which sug- 
gests that an axisymmetric Galaxy is incompatible with the MACHO results 
This must be a tentative conclusion, however, until one has taken into account 
the effect on a body's circular-speed curve Vc{r) of the body being strongly flat- 
ted. Binney, Bissantz &: Gerhard (2000) show that it is possible to choose the 
Galaxy's radial density profile in such a way that the required optical depth is 
obtained without generating a value of Vc that conflicts with observation. How- 
ever, such radial density proflles require more matter in the solar neighbourhood 
than observations of the Oort limit (Creze et al, 1998; Holmberg & Flynn, 2000) 
and the mass within 1.1 kpc of the plane (Kuijken & Gilmore, 1991) imply. Con- 
sequently, we can safely conclude that an axisymmetric Galaxy cannot have as 
large an optical depth as that reported. 

Can one achieve a higher optical depth within a give mass budget by making 
the bands elliptical rather than circular? Imagine deforming an initially circular 
band into an elliptical shape while holding constant the radius r at which the 
line of sight to our sources cuts the band. It is straightforward to show that if 
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the column density through the band to the sources is to be independent of the 
band's eccentricity e, its mass M(e) must satisfy 

M(e) = M(0)';f:°;'^ (17) 

where <j) is the angle between the band's major axis and the Sun-centre line. For 
(j) < IT /A, M{e) is a minimum with respect to e at 



'2 -secret). (18) 

Substituting equation (p^) in equation ( |T7|) we find the minimum mass to be 

Mniin = M(0)sin2(/., (19) 

and to require axis ratio gmin = tan0 (Zhao & Mao, 1996). 

For (f) = 20°, a value favoured by Binney et al. (1997), gmin = 0.36 and 
Mmin/M(0) = 0.64; for (/> = 15°, we find Qmin = 0.27 and Mmm/M(0) = 0.50. 
Hence, making the bands elliptical realistically reduces the mass required to 
generate a given optical depth by at most 50%. In practice we cannot reduce our 
requirement for mass by so large a factor because the structure of the Galaxy's 
stellar bar is strongly constrained by both near-IR photometry (Blitz &: Spergel, 
1991; Bissantz et al., 1997) and radio- frequency observations of gas that flows in 
the Galactic plane (Englmaier k Gerhard, 1999; Fux, 1999). Binney et al. (2000) 
estimate the possible reduction in mass by assuming that material at i? < 4 kpc 
forms a bar of optimal eccentricity whose long axis makes an angle of 20° with 
the Sun-centre line, while the Galaxy is axisymmetric at i? > 4 kpc. They show 
that if the vertical structure of such a Galaxy is chosen to be optimal for lensing, 
then a radial density profile can be found that nowhere exceeds the observed 
value of Vc and is also compatible with the constraints on the density of matter 
at Rq. They argue, however, that such Galaxy models can be excluded for two 
main reasons. First, they predict values of Vc that are too small at i? 350 pc 
because at small R they place significant mass high above the plane, where it 
contributes to r but not Vc- Second, these models do not leave enough room 
in Vc for (i) departures from the optimal vertical profile and (ii) the presence of 
matter, such as interstellar gas and non-baryonic matter, that contributes to Vc 
but not r. 

Hence, even though we don't know much about the Galactic force-field, we 
know enough to exclude the measured optical depth in Baade's window! In fact, 
if we were to take seriously the prediction of simulations of the cosmological 
clustering of CDM that dark matter contributes substantially to the mass inte- 
rior to the Sun (e.g., Navarro & Steinmetz, 2000), our predicted optical depth 
for Baade's window would be significantly less than 10~^, a factor of 4 or more 
below the measured value. Something is seriously wrong here. 
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